Path Integrals over Measurement Amplitudes: 
Practical Quantum Foundations for Signal Processing and Control 
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It is shown that classical control diagrams can be mapped one-to-one onto quantum path inte- 
grals over measurement amplitudes. To show the practical utility of this method, exact closed-form 
expressions are derived for the control dynamics and quantum noise levels of a test mass observed 
by a Fabry-Perot interferometer. This formalism provides an efhcient yet rigorous method for ana- 
lyzing complex systems such as interferometric gravity wave detectors and magnètic resonance force 
microscopy (MRFM) experiments. Quantum limits are conjectured for the sensitivity of interfero- 
metric observation of test mass trajectòries. 
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I. INTRODUCTION 

As reviewed by Mensky jïj, the formalism of path inte- 
grals over measurement amplitudes was first suggested by 
Feynman ||] and was subsequently worked out in greater 
detail by Mensky |, f§ and by Caves § §. 

We show in this article that the formalism of path inte- 
grals over measurement amplitudes can provide practical 
quantum foundations for control theory, and we illustrate 
these foundations by the worked example of resonant in- 
terferometric observation of a test mass. 

From a physics point of view we will work everything 
backwards. We will start, rather than finish, with a block 
diagram that describes the dynamics of a classical sys- 
tem that is subject to closed-loop control. We will show 
that such diagrams can be mapped one-to-one onto path 
integrals over measurement amplitudes. Then we will 
illustrate the physical and control-theoretic significance 
of each path integral term by analyzing a test mass ob- 
served by resonant optical interferometry. Finally, we 
will suggest that the dynamical behavior of such systems 
is connected to unsolved problems in quantum signal pro- 
cessing and cryptography. 
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FIG. 1: The control theory block diagram associated with 
equations (p~a|— c) , or equivalently the path integral (fïa|-b) . 
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II. CONTROL THEORY FOUNDATIONS 

Engineers commonly formulate control theory in terms 
of block diagrams and signal- flow graphs |7j]. A block 
diagram is conceptually similar to a Feynman diagram: 
it is a graphical representation of a set of equations. 

Figure shows a block diagram for a test mass whose 
position q{t) is continuously measured and controlled. 
The control force f c (t) is determined from the com- 
manded position y c (t) and the estimated position y(t) 
via kernels T and V: 



f c (t)= I dt'[T(t-t')(y c (t')-y(t'))] 

dt'T'(t-t')y(t'). (la) 



Here by convention the feedback has been separated into 
a control kernel T that is causal and a backaction kernel 
r' that is anticausal. 

Anticausal backaction kernels are a natural idiom in 
control theory. To see why, consider a present fluctua- 
tion in y(t') that represents a photon detected at time 
t', having bounced off the test mass at past times t. It 
follows that the backaction kernel T'(t — t') must be an- 
ticausal, in order to describe the past-time force exerted 
by those bounces. There is no implication that the back- 
action physics is noncausal. 

To anticipate, all the other kernels of Fig. [ïjare explic- 
itly causal, and so is the path integral that we will derive 
for the system dynamics. The overall formalism is therc- 
fore explicitly causal, as we will discuss following (||). 

The controller kernel T is the main subject of control 
theory; it typically serves some useful purpose like mov- 
ing the test mass to a commanded position or stabiliz- 
ing the system dynamics. Such kernels can be designed 
for optimal performance and essentially any desired 
causal kernel can be realized by digital technology || . 

The estimated position y(t) is determined from the test 
mass position q(t) and the measurement noise q n {t) via 
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the measurement kernel £ : 



y(t) 



dt'at-t')(q(t') + q n (t')). (lb) 



It follows that q(t) cannot be observed directly, but 
rather must bc estimatcd from y(t); such estimation plays 
a central role in control theory. 

Finally, the dynamical behavior of q(t) is determined 
by the external force f c (t), the process noise / n (í), and 
the control force f c (t) via the process kernel G: 



dt' G(t - t')(f e (t') + Uf ) + fc(t')). (lc) 



Our main statistical assumption is that q n (t) and f n {t) 
are stationary zero-mean random processes; we will show 
that path integrals naturally generate quantum noise 
with this property. Without loss of generality, we can 
further specify that Çn(í) and fn{t) are statistically in- 
dependent; for any given block diagram this can be ar- 
ranged by "pulling" correlated noise through G and suit- 
ably redefining q n (t) and / n (í). 

Then the Fourier-domain solution to (|ïa|-c) is such 
that y(tij) has mean value 



ÉH<?(W) (/e(w) + f(w)ífc(w)) 



1 + Í{lü)G(üj)(T{lü) + T<{lj)) 
and spectral density 

IÇHP^H + IGHi^/H) 



Sy{w) 



\l+£{Lj)G(Lü)(f{Lü)+r{Lü))\ 



(2a) 



(2b) 



The form of this result reflects our convention that the 
measurement noise q n (t) and the process noise / n (í) are 
uncorrelated. Here our notation and normalization con- 
vention for Fourier transforms and spectral densities is 



(3a) 
(3b) 
(3c) 



Sf(u)= / dr(/(0)/(r))e— , 

J — OO 
p OC 

K<»)= rfr/(r)e-^, 



dràr' iJ)(t, t ) 



r\ — i(«r+wV) 



with (. . .} designating an ensemble average. Thus our 
spectral densities are "two-sided." We regard the kernels 
of ([ïa|-c) as defined for all times r with T(r) = £(t) = 
G(r) = for r < and T'(t) = for r > 0. 

As a mathematical point, the functions / n (í) and ç n (í) 
in the block diagram of Fig. [ï] do not appear as indepen- 
dent functions in (|2a|— b) . Neither will these functions ap- 
pear as independent variables in our path integrals — not 
even as dummy variables of integration — nor will they 
appear in any subsequent part of our formalism. Their 
sole role is as mnemonic aids: they remind us to include 



Since f n (t) and q n (t) do not appear in our formal- 
ism as indendent functions, we have no mathematical 
or physical basis for assigning independent meanings 
to them. We can only speak of them in terms of a 
unitary equivalent noise having, e.g., spectral density 
S f (u>) + \G(Lü)\- 2 S q {ui) at the G block input, per Fig. |. 
Again anticipating future results, our sole motivation for 
maintaining Sf(uj) and S q (uj) as separate densities is to 
express the noise reciprocity relation (^9|) in a device- 
independent form. 

This unitary point of view is consonant with informa- 
tion theory, since in light of the above discussion it is 
not possible — even in principle — to infer independent vàl- 
ues for / n (í) and q n (t) from the measured quantity y(t). 
Furthermore, maintaining a unitary point of view will 
forestall conceptual difhculties in Section VI, where we 



compare path integral results with analyses of shot noise 
and radiation-pressure noise in the literature. 



III. PATH INTEGRAL REPRESENTATIONS 
OF CONTROL THEORY 

Now we seek a path integral that reproduces ((pa[-b)). 
We approach this as a purely mathematical exercise 
whose sole requirements are tractability and generality. 
Adopting the path integral notation of Brown [jïo]], we 
consider functionals of the form 



P(y(t)\fe(t),y c (t)) = 

[dq]exp [A(q(t),y{t)J c (t),y c {t))] 



(4a) 



Here P(y(t)\f e (t),y c (t)) is a Gaussian probability func- 
tional whose mean and variance must reproduce (|2a|— b) . 

For P to be of the required Gaussian form, the action 
functional A must be biquadratic in q(t) and y(t) and 
bilinear in f e (t) and y c (t). Note that q(t) appears only 
as a dummy variable of integration that is not yet iden- 
tified as the test mass trajectory of Fig. [ï| Adopting a 
conventional form that facilitates subsequent connection 
to control theory, the most general functional A that re- 
produces (pa]-b) can be written as 



A = 



Tíf + Tir 
Mv + Me - 

ÍMç] 



Lagrangian action 



external force and control 



backaction effects 



measurement amplitude (4b) 



Sf and S q in (2b). 



where {£, Tíf, Tir, -Mr 1 , Mg, M^, A4ç} are real-valued 
functionals that are given explicitly in Table [| in terms 
of measurement kernels {V , T, 6, ip, £}. Carrying through 
the path integral by methods that are essentially alge- 
braic JÏÏJ, we connect the kernels of Jibl ) to the control 
dynamics of (pa[-b), as summarized in Table pï|. 

This completes our goal of establishing a path integral 
representation of the control diagram of Fig. 0. 
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TABLE I: Table of functional densities that appear in (fïa|). See Table ^ for the connection between the kernels {r, V ,8, ip, £} 
of these densities and the system dynamics specified by Eqs. (pa[-b). 



Lagrangian test mass action 



External force and control functionals 

H f (q(t),f e (t)) = -q(t)f e (t) 



Hv{q{t),y{t),y c {t)) = -q{t) dt' T(t - t')(y c (t') - y(t')) 

•/ — oo 

Measurement functionals 

M T , (q(t),y(t)) = +q(t) H dt' V (t - t')y(t') 

J — oc 
[• oc 

M e (q(t), b(t)) = -b(t) / dt' 6(t - t')q{t') 

J — oo 



M^{q(i), r(t)) = +- / dt' dt" V(í - 1', t - t")q{t')q(t") <S> 



Mdq(t),y(t),r(t))-- 



ti í í°° \ 2 

.^{y(t)-j dt'i{t-t')q{t')) 



{The Lagrangian of a test particle of mass m with 
(optional) spring constant mtJo. 



{Generates the dynamical effects of the external 
force / c (í). 

{Generates the control force appropriate to the 
commanded position y c (t). 

{Generates a fluctuating backaction force that is 
deterministically correlated with fluctuations in y(t). 

{Generates a backaction force parameterized by the 
deterministic function b(t). 

{Generates a frequency-dependent backaction 
spring constant. 

{Correlates the measured value y(t) with the test 
mass position q(t). 



TABLE II: The connection between the measurement kernels 
{i/>, 8, Ç} of Table and the system dynamics and noise spec- 
ified by Eqs. (|a]-b) and Fig. |ï} 

Backaction dynamical effects 

G _1 (oj) = m{uj 2 — lo 2 ) + ipi 1 - , — w) <^ {dynamical kernel. 

f e {uj) = b(u>)8(— uj) <^> {backaction force. 

Measurement noise PSDs 



S q {u)-- 
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líMI 2 
h 2 



{measurement noise. 
& {force noise. 



IV. THE MEASUREMENT AMPLITUDE 
FOR OPTICAL INTERFEROMETRY 



To show what is gained by attacking the problem in this 
systematic way, we now calculate the measurement am- 
plitude kernels {V , V, 0, ip, £} for a resonant optical inter- 
ferometer. Then we systematically read off the system 
dynamics and quantum noise from Tables g and iïfl and 
Eqs. (|2a|-b). 

To carry through this calculation — and indeed to 
carry through any path integral/measurement amplitude 
calculation — it suffices to specify the classical optical 
scattering amplitude and the photon detection statistics. 
From classical physics we know that for a general single- 



port optical device the outgoing amplitude a ou t at time 
t is causally conditioned upon the past history of the in- 
ternal coordinate q(t) as follows: 



Q-out 



(t\q(t)) = a in e lÇ (l + í dt' a(t - t')q(t') 

\ J — OO 

+ / dt' [ dt" I3(t-t' 1 t~t")q{t')q(t")). (5) 



oo J — oo 



Here a- m is the input light amplitude (assumed constant), 
a and (3 are scattering kernels, and e 1 ^ is an overall phasc. 
Perturbations of order q 3 and highcr are neglected, and 
by convention we normalize a- m such that the photon in- 



put rate is r- u 



Causal boundary conditions are 



imposed. Then the mean rate (r(t)) at which photons are 
detected at time t is a functional of the past trajectory 
q(t): (r(t)) = |a out (í|ç(í))| 2 . 

Now we are ready for the key element of our formal- 
ism. We introduce as an ansatz the following fundamen- 
tal relation between the quantum measurement action of 
(pïa|--b) and the classical scattering amplitude (||): 



exp 



M(r(t),q(t)) 



<w(í|<7(í)) 



r(t) 



x exp 



a ut(t\q(t)) 
(r(t)-\a ont (t\q(t))\^ 



Aa s \a out (t\q(t))\ 



(6) 



Here M(r(t),q(t)) = M r < + M e + M,p + iM í is the 
measurement action that appears in (fïa]-b). 

We will present no ficld-theoretic justification for this 



ansatz, and in Section VIII we will present reasons for 
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TABLE III: Rules connecting the photon sideband amplitudes 
{a, ff\ to the measurement amplitudes {tp, 8, Ç, T} of Tables | 
and [o] for the interferometer configuration of Fig. ^| 



I- ç(í) 



Fields and variables 

y(u>) = r(u>) — 2-Kri n S(iü) 
Í>(üj) = 2nri n S(ui) 

-y 2 =a s r in <^> 
Measurement kernels 



í photon flux 
Ifluctuations. 

{light pressure. 

{Sets the flux spectral 
density S r ~ a a n n . 



^ r~/ \ -*/ m * íbackaction control 

fcouples light pressure 



Ikernel. 



6»(w) = -r'(-w) 

^(w) = — Sri n Im[2^(w, — w) -*4> 



Ito the test mass. 

{frequency-dependent 
spring constant. 

=rt n [a(w) + à*(— w)] {measurement kernel. 



thinking that a rigorous field-theoretic justification would 
involve deep quantum-informatic issues. Instead, our 
limited goal in this article will be to show that the ansatz 
reproduces, within a path integral formalism, known clas- 
sical and quantum physics. 

From the ansatz we immcdiatcly obtain rules for trans- 
lating the optical kernels a and j3 into the measurement 
amplitude M. of @. No physical insight is employed; 
instead we impose the purely algebraic requirement that 
the rules convert (JsJ) into an identity (up to 0{q 2 ) in 
the path integral functional). The rules are summarized 
in Table |ï|. Because the optical kernels a and /? sat- 
isfy causal boundary conditions, the measurement ampli- 
tude (||) and the resulting path integral (|4a|) are explicitly 
causal as promised in the discussion following (|la|). 

Non-rigorously, the ansatz can be derived by con- 
structing a numerical wave function simulation along the 
lines given by Gardiner and Zoller |TÏ| , with each pho- 
ton detected separately and accounted numerically. Such 
simulations become exponentially slower as the number 
of photons and the múltiple reflections of each photon are 
increased; this illustrates the well-known nonpolynomial 
difficulty of quantum simulation in general. Coding such 
simulations with a view to making thcm numerically effi- 
cient leads naturally to a path integral formalism. This is 
the path the author followed to the results of this article. 

The ansatz embodies two key physical principies, 
which arc both well satisfied in optical interferometry. 
First, individual photons are assumed to be detected dis- 
cretely, such that the decoherence associated with each 
detection event procedes to completion within a very 
short time compared to all other dynamical time scales 
of the system. Under these circumstances a well-defmed 
phase can be associated with photon detection. The in- 
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FIG. 2: A Fabry-Perot cavity, with the end mirror serving as 
a test mass with coordinate q(t). 



terference of these phases can be readily observed, which 
is of course the reason such measurements are called "in- 
terferometric" . The ansatz functional (a ou t/|a out |) r( ' í ' ) 
accumulates the test mass phase from repeated photon 
detections; this phase creates the quantum backaction. 

Second, it is assumed that large numbers of photons 
are detected, in which case the detection statistics can 
reasonably be described by a counting formula of the 
usual Gaussian form (the lower right-hand term in (^|)), 
such that the photon flux spectral density is S r = cv s ri n , 
with due allowance for photon number squeezing as 
parametrized by a s . The physical role of the Gaussian 
term is to restrict the domain of path integration, condi- 
tioned upon the flux measurement r(t), in precisely the 
manner envisioned by Feynman, Mensky, and Caves. 



V. A WORKED EXAMPLE: 
SINGLE-PORT FABRY-PEROT 
INTERFEROMTRY 

To make the path integral/measurement amplitude for- 
malism come alive we will apply it to an engineering anal- 
ysis of the Fabry-Perot cavity shown in Fig. g. This is a 
single-arm interferometer with single-port detection; it is 
not intended to represent a realistic gravity wave detec- 
tor. However, even this simple design exhibits complex 
dynamical and noise phenomena; our goal is to show how 
to use path integral methods in analyzing this behavior. 

Many of the results that we will obtain by path in- 
tegration have also recently been obtained by operator 
methods in the literature on gravity wave detection; this 
literature is reviewed in Section [VÏ]. We will find no sc- 
rious conflict between operator methods and the path 
integral/measurement amplitude formalism. 

In applying the path integral formalism, our sole com- 
putational job is to calculate the optical kernels à(uj) and 
/3(cj, — lu) of (^) for the Fabry-Perot cavity of Fig. B; the 
rest is substitution into the rules of Tables || and jiïj. 

The input light is right-going, as shown in Fig. H, and 
our phase convention is that it has space-time depen- 
dence a(x,t) = ai n e ïfc;r ~" ot , where the wave number k 
and the optical carrier freqency ljq are positive quantitics. 
This accords with the quantum mechanics phase conven- 
tion that a right-going quantum carrying momentum hk 
has a spatial wave function oc e íkx . With this convention, 
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positive mirror displacements afí) correspond to longer 
cavity lengths, as shown in Fig. |2], and positive forces act 
to push the mirrors apart. 

This same phase convention, when conjoined with the 
Fourier convention (3b), prescribes that optical side- 
bands at a frequency u m have a Doppler-shifted time- 



dependence e 



Positive-frequency sidebands 



(u m > 0) are therefore associated with redshifted optical 
quanta. This unintuitive convention will be important 
later on when we check energy conservation. 

Inspection of Tablc III shows that the dynamical and 



noise behavior of the system is completely determincd 
by the Fabry-Perot sideband amplitude à(u) and car- 
rier amplitude /3(w, — ui). A straightforward perturbative 
calculation yields for the sideband 



à{uj) 



2ike i(3u 



1 sin 2 p 



( e 2i(W-0) _ C0S p)(l - e-^cosp) ' 



(7a) 



and for the carrier 



à(u>)à(—u>) 



sin 2 p ■ 



2i sin(2ç!>) cosp 



2 sin p 

Here the power reflectivity of the input mirror is by def- 
inition cos 2 p, and the one-way optical phase length cf) of 
the cavity is cj> = kL + tt, with L the cavity length. The 
extra tt in the definition of <f> is conventional; it ensures 
that tuning to 4> = (mod 2ir) yields maximal intra- 
cavity optical power. For fixed 0, and therefore fixed 
optical power, we can maximize the sideband amplitude 
by tuning the mirror modulation frequency uj to ojt = </> 
(mod 27r) ; this is physically equivalent to tuning the side- 
band on-resonance. 

The same perturbative calculation yields for the phase 
of the output light 



e 2l ^ cos p 



1 



-2i<p 



(8) 



cos p 

and it is easy to check that |e^| = 1. It is also useful to 
know the cavity power gain Ç ((/>): 

2 



Q{<P) 



d .f cavity power 
input power 



sm p 



cos p 



-2i<j> 



(9) 



Q(4>) is a period-7r function of <j), peaked about 0^0, 
and having half-width-half-maximum (f>yr 



- \J cos p + cos 1 p — 2. 



(10) 



By construction, <pjr is the extra phase length required to 
reduce the on-resonance cavity power by 1/2. Conven- 
tionally T = tt / {2<j)jr) is called the finesse of the cavity, 
and for high-finesse cavities the on-resonance intracav- 
ity power gain is Ç7(0) ~ 2T jix = <t>^r l ■ We will always 
specify changes in cavity length as múltiples of <j>p. 

Now we are ready to specify an interferometer design. 
We choose length scales and power levels that are char- 
acteristic of recent proposals for advanced gravity-wave 
detectors 12 : 




-15 -10 -5 5 10 15 

cavity length (<j>/<pjr) 



-15 -10 -5 5 10 15 

cavity length (<p/(pjr) 



FIG. 3: Static force and spring constant generated by the 
intra-cavity light. Here <j> = f kL + n, with L the cavity length 
and k the optical wavenumber, and <\>t is defined in (flOh. 



cavity length: L 
detected power: P out 
on-resonance power: P max 
light wavelength: 

test mass: m 



= 4 km 
= 180 W 
= 830 kW 
A = 1064 nm 
40 kg 



(7b) This is a high-finesse cavity, with T ~ 7240, correspond- 



ing to a mirror power reflectivity cos p ~ 0.99913. 

We begin by considering the static behavior of the sys- 
tem. The path integral prediction for the light force / c (t) 
on the mirror can be read off from Tables O and [II : 



/e(í) 



00 du pi^t 



AH 



» duj 

-CK5 2-7T 

r°° du 



= -nnf'(0) 

•a*(o)] 



by definition 
by Table g 
by Table |TJ 
by evaluation 



fi r , , 

g [a(0) 



by Table III 



By explicit calculation we check that this accords with 
the Standard expression for light pressure: 



/ e (í) = 2hkr in Ç(cj ) ). 



(11) 



Here 2hk is the momentum transferred to the mirror by 
the reflection of a single photon and r- ln Q{(j)) is the flux 
of photons incident on the mirror. 

Similarly, from Table |ï| we see that the zero-frequency 
spring constant exerted by the light is fc S pring(0) = 
lim^^o il>(u>, — uj), and by explicit calculation we check 
that this accords with the spring constant predicted from 
the cavity gain GÍ4 1 )'- 



s (0) 



-2hk 2 r-, 



dG(4>) 



(12) 



As a final check on the static physics, it is straight- 
forward to show that the large T limit of k spj: i ng ((j}) 
precisely accords with the Fabry-Perot spring constant 
calculated by Braginsky, Khalili, and Volikov |Ï3| from 
classical physics. We thus confirm that the path inte- 
gral/measurement amplitude method accurately repro- 
duces known results relating to static optical forces and 
springs in resonant cavities. 
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In subsequent calculations we will not show all the 
steps, but our results are always obtained by a similarly 
direct application of the measurement amplitude rules 
given in Tables || |ï|. These rules are readily processed 
by symbolic programs; this reduces the incidence of al- 
gebraic error. Analytic continuation is straightforward 
because the kernels are given in closed form; we will see 
that this simplifies stability analysis. 

Now we turn our attention to the practical challenge 
of tuning the interferometer for dynamical stability and 
good noisc pcrformance. The static optical force and 
spring constant are shown in Fig. |^. The optical forces 
are weak — a few millinewtons at most — but the spring 
constant can approach 100 N//xm, which is extraordinar- 
ily stiff. For an optical beam of nominal diameter 20 cm 
and length 4 km the stiffness is equivalent to a modulus 
~ 12.3 TPa, which is twelve times stiffer than an equiv- 
alent bar of diamond. 

This illustrates that light itsclf can scrve as a structural 
material, as was first recognized and explored for dcsign 
purposes by Braginsky, Gorodetsky and Khalili 
|Ï6| and subsequently by Buonanno and Chen JÏ7] , 

By inspection of Fig. || we see that static stability is 
possible if and only if the interferometer is tuned "long" 
(i.e., cavity phase length > 0), and we henceforth con- 
fine our attention to this range. In practice, tuning is 
achieved by applying a few mN of force to the mirror to 
press it against the optical spring; the magnitude of the 
static force determines the equilibrium cavity length and 
therefore the optical tuning. 

Static stability does not guarantee dynamic stability, 
and on physical grounds we expect Fabry-Perot cavities 
to be dynamically unstable. We reason as follows: if we 
push the mirror against the optical spring, the spring will 
push back, but only after a time lag while the intracavity 
intensity builds up. By oscillating the mirror, we can 
continuously extract energy from the system. 

Physically, the possibility of energy extraction indi- 
cates the presence of dynamical instability. We will 
now prové that such instabilities exist — for all cavity 
tunings — by calculating the transfer function of the sys- 
tem in closed analytic form. 

We begin by noting that the time-averaged flux of 
output photons from any linear lossless optical device 
must equal the flux of input photons. It is easy to check 
that the Fabry-Perot amplitudes (f7a|-b) satisfy this con- 
straint, which requires that {a, /?} satisfy 

\à(uj)\ 2 + |<S(-w)| 2 = -2 * Re[j§(«, -to)] , (13) 
[sideband photon flux] = [decreased carrier flux]. 

Physically, photons that disappear from the carrier must 
reappear in the sidebands. 

This does not guarantee energy conservation, since 
the outgoing sideband photons are Doppler shifted per 
the discussion preceding ([faj-b). To check energy con- 
servation we apply an external force / c (í) to the mir- 
ror, such that the mirror is driven at an amplitude 
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FIG. 4: Dynamical behavior of the test mass. (a) Poles of 
the transfer function T(lú) as the cavity tuning is varied over 
< 4> < oo. The specific tunings (p/fa € {0.01,0.1,1,10} 
are marked with a '+'. A fixed zero at s = is marked 
with a 'o'. (b-c) A Bode plot of T(u) evaluated for the 
particular cavity tuning <j> — 10 (fr^. The magnitude axis is 
20 log 10 \f (277/) N/ M W|. 



q(t) — qo cos(w m í) . Then mechanical energy is supplied 
to the mirror at a rate given by (|2a| ) as 

(fe{t)q(t)) t = ÍWzn J [f '(-W m )|(-0, m ) 

-f'(u m )ï(iü m )\ (14) 

with (. . .) t denoting a time average. After substitution 
from Tables || and III this is 



(fe(t)q(t))t 



\à(-u m )\ 2 - \à(üj m )\ 2 ] . (15) 



Taking into account the Doppler sign convention dis- 
cussed at the start of this section, we recognize this as 
precisely the excess optical power emitted in the side- 
bands. Thus the optomechanical instability is energeti- 
cally driven by Doppler shifts in the sidebands, such that 
energy is explicitly conserved overall. 

Now we analyze the instability in detail. From (p^), 
the transfer function T(lu) for the cavity is 



f (w) = hkc 



hkc ■ 



i + e(w)G(w)r'(w) 



(16) 



where the factor hkc normalizes T(lü) to units of watts 
of detected optical power per newton of applied force. 

In control theory — where Laplace transforms are morc 
common than Fourier transforms — it is Standard prac- 
tice to plot the dominant poles and zeros of the transfer 
function in the complex s-plane, with s = iu). 

The path integral/measurement amplitude formalism 
gives T(u>) in analytic form, and when such forms are 
available, a Standard technique in control theory is to 
calculate a Padé approximant to T(u>) using a symbol 
manipulation program. Carrying through this calcula- 
tion, we find that a {1,4} approximant yields a good 
fit, and is comprised by a fixed zero at s = and four 
dominant poles whose trajectòries are shown in Fig. H. 
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The fixed zero at s = means that the cavity has zero 
static sensitivity at all tunings. Physically, this means 
that at zero frequency the measured output photon flux 
equals the input flux, no matter what the cavity length, 
as enforced by photon conservation (p"3|). 

Dynamically, the cavity is optomechanically unstable 
at all tunings, with the strongest instability at 4>/4>t ~ 1. 
To stabilize the cavity, an additional control kernel T(lj) 
must be added. In principle, a perfectly linear and noise- 
less control kernel will not altcr the signal-to-noise ra- 
tio [|[ ||] — because control kernels affect signal and noise 
equally — but nonetheless it is good engineering practice 
to choose a cavity tuning such that the control challenges 
are not too great. 

We choose a far-ofï-resonance cavity tuning 0tunc = 
1O0jf. Per (H), this tuning reduces the cavity optical 
power from the peak power of 830 kW to 8.2 kW — 
a 99% reduction. In the ordinary course of events, we 
might expect such low power to greatly diminish the sen- 
sitivity of the interfcrometer. 

However, the design compensates for the low cavity 
power by exploiting two mitigating factors. First, on 
physical grounds we expect that the signal sideband will 
be resonant with the cavity, and therefore passively am- 
plified, at a frequency w pticai determined by (f7aj) to be 
^o P ticai/(27r) ~ 0tuno/(27rr) ~ 25.9 Hz. Second, we ex- 
pect the system to exhibit a mechanical resonance at a 
frequency cj mech such that 4>(u> mec h, -c^Wh) - mu) L·ch^ 
i.e., at a frequency determined by the strength of the op- 
tical spring, and we further expect this mechanical res- 
onance to be unstable. For 0tune = 10 4>j^ a numerical 
analysis predicts this resonance at w mec h/(27r) ~ 23.7 Hz. 

Thus, on physical grounds we expect the transfer func- 
tion to exhibit one stable optical pole and one unstable 
mechanical pole, both at about 20 Hz. 

These expectations are in excellent accord with the 
Padé analysis shown in Fig. [|, which at <fi = 1O0jf ex- 
hibits a stable pole at 21.2 Hz with a quality Q = 1.90 
and an unstable pole at 15.9 Hz with Q = —2.59. These 
resonant freqüències are slightly reduced relative to the 
above rule-of-thumb expectations; this presumably re- 
flects damping effects (which characteristically lower res- 
onant freqüències) combined with the optomechanical 
coupling generated by the backaction kernel f". 

Bearing in mind the consensus view of control theorists 
that "generally speaking, an unstable system is consid- 
ered to be useless" , these results provide a well-posed 
starting point for addressing important practical qües- 
tions such as: is this Fabry-Perot system "observable" 
and "controllable" in the rigorous sense that these terms 
have in control theory? And if so, what would be a suit- 
able design for the control kernel r(w)? 

We note that in control theory — and in any continu- 
ous measurement theory — there is no sharp distinction 
between a "position meter" and a "velocity meter." If 
position is observable, then so is velocity, via a differen- 
tiating filter. Conversely, if velocity is observable, then 
so is position, via an integrating filter. 



We will not consider these control issues further be- 
cause an article-length exposition would be required, and 
because they are a Standard tòpic in control engineering 
textbooks [[7| . Instead, we will simply assume that a sta- 
bilizing controller is present, we will further assume that 
it contributes negligible noise, and we will proceed to an- 
alyze the sensitivity of the interferometer. 

In keeping with accepted practice of the gravity-wave 
community, we conflate all noise sources into single equiv- 
alent force noise having spectral density Sj ot (u>), and we 
express this net force noise as an equivalent strain noise 
5£ ot (u;) according to 



Çftot 



S} ot (w) 
{muj 2 L) 2 



(17) 



where L is the arm length and m is the test mass. 
Physically, this convention acknowledges that audio- 
frequency gravity waves, when observed over kilometer 
length scales, are dynamically equivalent to tidal forces. 



Combining this convention with (2b), we obtain for the 
total equivalent strain noise 



Cftot 



S q (u) 



g/M 

{mu 2 L) 2 (muj 2 L) 2 \G{üj)\ 



(18) 



process 


+ 


measurement 


noise 




noise 



where the functional forms of Sf(ui), S q (ui), and G(u>) are 
given in Tables and 0. The feedback kernels {f ', f, £} 
do not enter because they affect signal and noise equally. 

Even though q n (t) and f n (t) are statistically indepen- 
dent, the above result is fully consonant with predictions 
from field theory (jÏ7| |Ï8|, |ï^, [2(| of correlations between 
"shot noise" and "radiation-pressure noise." We identify 
shot noise with fluctuations in y n (t) and radiation noise 
with the sum f n (t) + / c (*)< Then from Fig. [ï] we see that 
the combined effect of the measurement kernel £ and the 
backaction kernel V is to apply a fluctuating radiation 
force f c (t) that is deterministically related to the shot 
noise y n (t). This is the path integral mechanism that 
correlates shot noise with radiation-pressure noise. Note 
that /c(í) does not enter into the total strain noise ( |Ï8| ) 
because it is already accounted by the spectral density 
S q (uj) via the feedback kernels £ and V in Fig. 

Figure shows the resulting performance of the ref- 
erence design. As predicted by the Bode diagram of 
Fig. |], the sensitivity of the system is greatly amplificd 
in the frequency band where mechanical and optical res- 
onances coincide, viz., the band 10-30 Hz. From Table 
we find that the process noise Sf(ui) and measurement 
noise S q (u>) satisfy a device-independent equality 



S f (u)S q (w) = h 2 /A. 



(19) 



Minimizing 5jj ot (w) subject to this equality defines the 
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FIG. 5: Noise performance of the Fabry-Perot cavity. The 
rms total strain noise is defined to be [2S'j; ot /(l second)] 1//2 , 
and all the other noise densities are normalized similarly. The 
equivalent one-sided bandwidth is one hertz. 



sprung mass quantum limit S^ pruns (oj): 



sr uns M=\™* «TM 

{Sf,Sq} 



i 2 uj 4 L 2 \G(uj) 



(20) 



where the dynamical kernel G(u>) is held fixed during 
minimization. It follows that S^ pTUTls (uj) sets a rigorous 



lower bound to the interferometer noise S^ ot (ui) : 



(21) 



This lower bound is included in Fig. || and is seen to 
be saturated at two discrete freqüències: ~ 15 Hz and 
~ 25 Hz. Physically speaking, at these freqüències the 
process noise Sf and the measurement noise S q are opti- 
mally balanced for strain detection. 

As a final check, the assumption of free test mass dy- 
namics yields the Standard quantum limit (SQL) 



jSQL, 



lim 

1/|G|-»TO 

h 

muj 2 L 2 



Çftot 



(22) 



in accord with the literature 1 . 

As shown in Fig. 0, the reference design beats the SQL 
in the 10-30 Hz band. But this does not signify any 
evasion of the rigorous quantum limits (pj|~pl·|), because 



This agrees 



1 The literature value is St^ L (a>) = 8h/ {rn. m i t ü) s 
with (E2J) upon taking m m ; r — > 4m mo( j c , with m mo[ j e the mo- 
tional mass of a four-mirror interferometer, then m moc [ c — » m, 
with m our single-mirror mass, and finally inserting a factor of 
1/2 to convert to a two-sided spectral density. 



TABLE IV: Variables in which quantum fluctuations are 
commonly identified. 



Íin(í) 


the làser source current, 


A in (í) 


the gauge fïeld at the input port(s), 


■A C a.v (i ) 


the gauge fïeld within the cavity, 


jmir(í) 


the mirror current, 


Aout(í) 


the gauge field at the output port(s), 


jout(í) 


the photodiode sink. 



the SQL assumption of free test mass dynamics is not 
justihed, viz., the dynamical kernel G(u>) differs greatly 
from the free kernel — l/(mw 2 ) in consequence of optical 
forces and springs. 

These results illustrate a fundamental principio of field 
theory: all physics can be derived from the scattering 
matrix. In our case the scattering matrix is the optical 
amplitude (|7a|-|7b"|) , and from this sole input the path inte- 
gral/measurement amplitude formalism constructs both 
the classical dynamics and the quantum noise. 



VI. ACCORD WITH THE GRAVITY WAVE 
DETECTION LITERATURE 

Now we will show that the path integral results of the 
preceding section — both dynamical and noise-related — 
accord with prio r results from the gravity-wave detection 
community || || || ^ |lj || || || 

p5| , which were obtaincd mainly by operator and frcld- 
theoretic methods. 

Showing accord is daunting because the gravity-wave 
community has a tradition — extending back at least 
twenty years — of "lively controvèrsies" pï[ . 

Much recent discussion has been stimulated by the 
work of Braginsky and colleagues PH , who have been 
cited |0 as showing that "the test-mass wave-function 
aspect of the uncertainty principle is irrelevant to the 
operation of a [gravity-wave] interferometer" . 

To put a sharp point on the issue, how can the path 
integral/measurement amplitude formalism, in which the 
test mass is explicitly quantized but light is not, be equiv- 
alent to other — seemingly opposite — formalisms in which 
the test mass is not explicitly quantized but light is? 

In showing that there need be no contradiction, we will 
build upon a seminal article by Caves \2l and an analy- 
sis of photodetection by Gardiner and Zoller JÏJ]. Caves' 
analysis revealed "two different, but equivalent points of 
view regarding the origin of . . . radiation-pressure fluctu- 
ations." Seeking further equivalent points of view, we can 
idcntify in Gardiner and Zollcr's analysis (and in much 
other quantum òptics literature) at least six variables in 
which quantum fluctuations occur (Tablc IV). 

Viewed as field operators, these variables are linked by 
Maxwell's equations, such that fluctuations in any one 
operator determine the fluctuations of all the others (up 
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to boundary conditions on the optical field). Thus, any 
one of these variables can reasonably serve as the focus 
of an "equivalent point of view" in the sense of Caves. 

The path integral/measurement amplitude method 
amounts to a point of view that is focussed upon the test 
mass trajectory q(t) and can be formalized as follows: 

• The quantum dynamics of the mirror are embod- 
ied (pa|-b) in a path integral over q(t). 

• MaxwelPs equations are enforced by the optical 
kernel (|5|), with q(t) as the source term. 

• Optical boundary conditions are specified via the 
measurement amplitude (^|). 

Adhcrants of this point of view can noncthclcss consis- 
tently agree with the very different point of view of Caves' 
articles ^ï), |22|| , which focus upon vacuum fluctuations 
at the input ports. In a measurement amplitude formal- 
ism the port fluctuations appear implicitly in the photon 
counting statistics of the measurement am plitude (|), m 
accord with Gardiner and Zoller's dictum jO]: 

Under the conditions that normally apply for 
a practical photodetector, the 'out' electron 
field has the same statistics as the 'in' photon 
field. 

Thus, vacuum fluctuations entering at the input port nec- 
essarily appear in the photon statistics. This reconciles — 
at least in principle — the port-oriented point of view with 
the path integral/measurement amplitude formalism. 

This suggests the general principle that a path inte- 
gral/measurement amplitude analysis should agree with 
any other analysis that treats at least one variable 
quantum mechanically, enforces Maxwell's and Newton's 
equations, and imposes compatible boundary conditions 
on the optical fields. 

As a test of this "many viewpoints" principle, we have 
carried through a path integral analysis of each of the 
measurement schemes that were analyzed in [^5| by op- 
erator methods; we find exact agreement between the 
two formalisms in all cases. This work will be reported 
in a separate article. However, a conceptual issue arose 
in which the language of control theory proved more pre- 
cise than the language of physics; this precision played a 
key role in reconciling the two viewpoints. 

The issue is: what is a free mass? To a control engineer 
the question is ill-posed, because an appropriate control 
kernel can create dynamics that arc equivalent to a vir- 
tual spring, even though no physical spring is present. 
Conversely, a physical spring attached to a test mass can 
be veiled by a control kernel, such that the controlled 
dynamics are equivalent to those of a free mass. 

Veiled springs pose a conceptual obstacle in measure- 
ment theory because they allow violation of the free-mass 
Standard quantum limit by test masses that are only 
seemingly free. We found in |25| sever al examples of me- 
ters that — upon computing an equivalent measurement 
amplitude, path integral, and control diagram — proved 
quantum-mechanically equivalent to a physical spring 



plus a spring-veiling controller. In every case the reci- 
procity relation (19) was satisfied, such that violations of 
the Standard quantum limit were due to the veiled spring. 



VII. CONJECTURED LÍMITS TO 
INTERFEROMETRIC TEST MASS 
OBSERVATION 

We set forth in this section conjectured limits that, if cor- 
rect, constrain all designs for interferometric gravity wave 
detectors, including past and future quantum nondemo- 
lition designs. To maintain an intel·lectual equilibrium — 
and to help sustain the gravity wave community's tra- 
ditio n of "lively controversy" — we will suggest in Sec- 
tion 



VIII several approaches by which these conjectures 



might be proved wrong or evaded. 

To start, we propose the following conventional dcfini- 
tion of a free mass: 

Definition 1 A free mass is defined to have a transfer 
function G{uj) — q(u))/ f(u>) — —l/(muj 2 ) = l/(ms 2 ). 

Here s is the Laplace variable traditionally preferred over 
lo by control engineers, and Q(u>) includes both mechan- 
ical and optical springs. This definition unveils hidden 
springs, and from a control engineering point of view is 
the most natural definition. 

We then propose the followingconjecture, which for- 
malizes the path integral result (|19|): 

Conjecture 1 For any stationary lossless interferomet- 
ric measurement processes, the measurement noise spec- 
tral density S q (uj) and the process noise spectral density 
Sf(u>) satisfy an exact equality 



S q (u,)S f (u) = h 2 /A 
and these noise processes are uncorrelated. 



(23) 



Here the "stationary" constraint excludes stroboscopic 
]Ï5[ and squeezed j2í| measurements, and "lossless" is 
understood to mean "no unobserved decoherence." We 
recali from the discussions following (^tj) and (|ï^) that 
the measurement noise q n (t) and process noise / n (í) are 
not equivalent to shot noise and radiation-pressure noise, 
but rather are mnemonic aids whose sole role is to remind 
us to include Sf and S q in (2b|). 

The point of Conjecture [ï is to suggest that ( p3| ) need 
not be regarded as a limit to be approached, but instead 
provides us with an exact law of nature that even the 
most clumsily designed experiments cannot violate, pro- 
vided only that all decoherence is observed and no ob- 
servations are discarded. This viewpoint facili tates the 



informàtic investigations we propose in Section VIII 



Definition [I] and Conjecture |ï| lead immediately to a 
lemma (derived in (|Ï£ - 20 )) : 
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Lemma 1 For any stationary interferometric test mass 
measurement, the spectral density S^ oí (w) of the equiva- 
lent gravitational strain noise satisfies an inequality 



múltiple output ports. A natural n-port ansatz is 



Sl°\uj) > 



h 



m 2 uAL 2 |GH| 



(24) 



where m is the reduced mirror mass, L is the arm length, 
lo is the angular observation frequency, and G(ui) is the 
test mass transfer function, including optical forces. 

Since this inequality — the sprung mass quantum limit 
holds even for squeezed photon detection statistics (as 
discussed following (^)), the point of Lemma |ï| is to sug- 
gest the strong hypothesis that all stationary measure- 
ment schemes for exceeding the Standard quantum limit, 
if analyzed from the path integral/measurement ampli- 
tude point of view, and with care taken to unveil hidden 
springs, are equivalent to the design strategy of the pre- 
vious section, which can be formalizcd as follows: 

• Install an optical or mechanical spring that in- 
creases |G(w)| relative to the free mass value of 
l/(mw 2 ), 

• Simultaneously tune the sideband response of the 
cavity to balance S a {ui) and Sh{<+>) such that the 
sprung mass limit (E4T) is saturated over the broad- 
est feasible bandwidth, and 

• Install a control kernel T(uj) to quench any op- 
tomechanical instabilities and — if desired — alter 
or veil the dynamical effects of the optical spring. 

The nondemolition meters proposed in p5| are consistent 
with this strategy, but there are many other proposed me- 
ters in the literature that remain to be considered before 
it could be considered general. 

If Lemma is correct, then optimizing the sensitivity of 
interferometric gravity wave detectors is a problem that 
can be posed purely in terms of classical optomechanical 
design. Because the theoretical and practical limits to 
maximizing G(u>) are not known, Lemma imposes no 
fundamental limit — quantum or classical — on the sensi- 
tivity of interferometric gravity wave detection. 



VIII. DISCUSSION 

Conjecture [ï] and Lemma |ï| are suggested by the path in- 
tegral/measurement amplitude formalism, but they are 
far from proved. We will now outline a program by which 
they might be proved wrong or evaded. Beyond its intrin- 
sic scientific interest, this program would advance at least 
three practical goals: quantum cryptography, single-spin 
imaging, and interferometric gravity-wave detection. 

Quantum cryptography and quantum entanglement 
considerations arise naturally when we consider how to 
generalize the measurement amplitude (^|) to the case of 



exp 



M({n(t)},Q(t)) 
ih 



■ni 

iE~L,n 



Oi(t\q(t)) 



\ai(t\q(t))\ 



n(t) 



x exp 



-( n {t) - \ ai (t\q(t))\ 



2 i 2 



4o i |o < (t|ï(i))|í 



(25) 



where {ai(t\q(t)),ri(t),ai} specify the amplitude func- 
tional, detection rate, and photon count squeezing at 
the i'th output port. This measurement amplitude — 
or a similar expression — would have to be rigorously 
grounded in field theory before Conjecture [ï] could be 
regarded as a theorem. Furthermore, the design analy- 
sis of real-world gravity-wave detectors, whether analyt- 
ically or by quantum numerical simulation, also rcquircs 
an explicit rt-port measurement amplitude. 

The following thought experiment suggests how chal- 
lenging such a field-theoretic grounding might be. Con- 
sider a four-port interferometer observing a single test 
mass, in which Alice monitors Ports 1 and 2 while Bob 
monitors Ports 3 and 4. 

Alice and Bob decide — independently and secretly — 
how to process their ports. For example, Alice can decide 
to count photon rates |eti| 2 and |a.2 1 2 , or alternatively she 
can count \ai — ia 2 \ 2 /2 and \ai + ia 2 \ 2 /2; Alice's data 
records will in general be very different depending on her 
choice, as will her inferred vàlues of S q (uj) and Sf(uj). 

Alice's choices must be invisible to Bob, and Bob's 
choices must be invisible to Alice; otherwise causality is 
violated. But depending on the quantum dynamics of 
the test mass, there is at least the possibility of quan- 
tum entanglement of Alice and Bob's port amplitudes. 
Furthermore, Alice and Bob have the option — at least 
in principle — of storing their light away, for analysis at 
some future time by a method to be decided later; such 
delayed choices must also be mutually transparent. 

Such thought experiments suggest that rigorously jus- 
tifying or refuting (p3|-p5|) will encompass nontrivial is- 
sues of quantum entanglement, consistent with a recent 
proposal by Marshall et al. [g7J . 

Quantum entanglement issues appear with redoubled 
subtlety when we consider magnètic resonance force mi- 
croscopy (MRFM). As with interferometric gravity wave 
detection, MRFM experiments monitor test masses by 
optical interferometry [ p8| p^ , |30| . Their diffcring phys- 
ical scale — nanograms, nanowatts, and nanometers for 
MRFM interferometers versus kilograms, kilowatts, and 
kilometers for gravity wave interferometers — is not par- 
ticularly relevant to the physics. More fundamentally 
different is the MRFM community's goal of observing 
the non-classical force signal from an individual spin. 

Early work in MRFM included the test mass quan- 
tum dynamics, but did not include a quantum analysis 
of the measurement process . Conversely, direct inter- 
ferometric observation of a single spin has been analyzed 
pï| , but this "toy" analysis did not include any test mass 
dynamics. Thus, no integral quantum analysis of a com- 
bined interferometer/spin/test-mass system is available 
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at present. As MRFM technology approaches attonew- 
ton force sensitivity |32| , such that long-envisioned single- 
spin detection and bioimaging applications |3(J approach 
feasibility, this fundamental quantum measurement chal- 
lenge is gaining in urgency. 

The statistical nature of the transition between spin-up 
and spin-down signals has crucial practical significance 
for the MRFM community; it strongly conditions the de- 
sign of optimal signal processing algorithms. This is is 
a practical embodiment of a decades-old question: when 
and how does a quantum wave function collapse? 

Similarly gaining in urgency is the practical challenge 
of how best to tune and operate gravity wave interferome- 
ters. For want of theoretical guidance, the interferomctcr 
of Section [v] was tuned empirically. Had we sent the out- 
put photons to Alice — per the discussion above — along 
with a homodyne reference, Alice might have achieved 
much better sensitivity, viz., a substantially more opti- 
mal balance between S q (u>) and Sf(u>). 

Alice's secretly improved sensitivity has to be trans- 
parent to Bob's simultaneous observation. Thus, part of 
Alice and Bob's communication challenge is to agree on 
how best to establish a consensus test mass trajectory, 
and how best to distinguish shot noise from radiation- 
pressure noise in their combined data records. The re- 
sulting Alice-Bob dialog would cast new light on these 
contcntious issues — doubly so if they were sharing non- 
classical spin signals in an MRFM context. 
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